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Abstract 

Motivated by the possibility to use high intensity neutrino beams for neutrino- 
nucleon scattering experiments we analyze charged current induced exclusive meson 
production within the framework of generalized parton distributions. The cross sec- 
~*ri \ tion for hard exclusive D s production is estimated in this formalism to leading order 

/"■s ■ in QCD. The integrated cross section proves to be sizable. It is shown that the con- 

sidered process is well suited to provide novel information on the gluon structure of 
nucleons, which is contained in the generalized gluon parton distribution. 

PACS numbers: 13.15.+g, 12.15.Ji, 12.38.Bx 

CC) ■ The factorization theorem [|], states that up to power suppressed terms every contribu- 

tion to the amplitude for hard exclusive meson production can be written as a convolution 
of a generalized parton distribution (GPD), a distribution amplitude, and a hard part. 
This has recently been applied for the investigation of electroproduction of single light 
^ mesons [|3| and meson pairs ||. The fact that in the near future also high intensity neu- 

Y" 1 ' trino beams might be available for lepton-nucleon scattering experiments |5|, || motivates 

the present study extending the formalism to charged current induced processes. 
<-| ! We consider the example 



X 



up + N -+ /i + + N + D~ . (1) 

Similar processes have already been subject of experimental studies J/J. For the analysis 
presented here, however, Q 2 has to be large compared to —t and the squared masses of 
the involved particles. The leading order amplitude is given by the sum of three diagrams 
involving a gluon GPD (Fig. |l]a and diagrams obtained by an interchange of vertices) and 
two diagrams with a contribution of the (polarized and unpolarized) strange quark GPD 
(Fig. |l|b plus one diagram with changed order of vertices). It is important to note that 
the only dependence on quark distributions in the nucleon comes in via the strange sea. 
Therefore the considered process is a good probe for the gluon GPD, which dominates 
the amplitude. To leading order in the strong coupling constant a s and neglecting terms 
of the order rn\)J(Q 2 + mjj ) 2 where rri£> s is the mass of the D s meson the amplitude for 
the subprocess 

Wl (q) + N(p)^ D- (q') + N(p') (2) 

is given by 

T = (N(p'),D;(q')\J ch ^ d -e L \N(p)) 
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Figure 1: Two of the five contributing diagrams. 
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(3) 



where the skewedness parameter £ is related to the Bjorken variable by £ = £bj(Q 2 + 
m 2 Ds )/{2Q 2 — xbj(<5 2 + mfjj). The mass of the charm quark m c is understood as pole 
mass, in the following a value of m c = 1.5 GeV is used. For the GPD's Fa(r,^,t), 
Fa(r,£,t), F s (t, £, t), and F s (r, £,t) the notation of Ref. || is used and the distribution 
amplitude for the D~ meson is defined by 
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(4) 



The differential cross section for leptoproduction is given in terms of this amplitude by 



da 



Xb\ 



Q 2 



dx Bj dQ 2 dt 2(4 7 r) 3 sm 2 # VF Q 2 (Q 2 + M 2 / ) 2 V 2x Bi p-l 
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Em 2 , 



(5) 
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Figure 2: The differential cross section for exclusive D~ production as a function of xbj or 
Q 2 respectively. The dashed lines show the contribution stemming from the gluon GPD 
F(r,£,t) alone. The small contribution of the polarized gluon GPD F(r,£,t) multiplied 
with a factor 100 is plotted with dotted lines. 

where / is the neutrino momentum. 

For an numerical estimate of the cross section we model the distribution amplitude 
following M as 
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DZ 



N. 



D t 



z(l — z) exp 
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D Sz 2 



(6) 



taking for the parameter u the value u> = 1.38 GeV obtained in [10J as the best fit for the 
D meson. The normalization constant No. has to be chosen such that the sum rule 



d* $£,-(*) =f Ds 



(7) 



is satisfied, where we adopt for the decay constant /d s the value fo a = 270 MeV as an 



average of the results obtained so far in lattice calculations [11], see also [12] for an earlier 
review. 

The gluon and quark GPD's are parameterized as in Q] combining Radyushkin's model 
||, [13| with the parameterizations of the usual (forward) parton distributions of Refs. [p 



(MRS (A')) and [T§ (Gluon A (NLO)). For the t-dependence of the GPD's we adopt the 
factorized ansatz F(r,£,t) = F(r,£,,0)F 9 (t)/Fg(0) and use the parameterization of |L(J 
for the gluon form factor Fg. For the strong coupling constant the two loop result is taken 
with N f = 4 and A^ D = 250 MeV. 

Figure |2] shows the results obtained for the differential cross sections dcr/d^Bj and 
da/dQ 2 where t = (p — p') 2 has been integrated over the interval £ m i n = m 2 N ^ 2 /(l — £ 2 ) < 
— t < 2 GeV 2 and the neutrino energy has been chosen as E D = 34 GeV. For the plot of 
the ^Bj-dependence Q 2 has been integrated from 12 GeV 2 to the upper bound given by 
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Figure 3: The differential cross sections from Fig. ^| compared with the results obtained 
for the asymptotic form of the distribution amplitude & D - (dashed lines) and by modeling 
the GPD's by their forward limit (dotted lines). 
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Figure 4: The result without mass corrections of the order m 2 D JQ 2 . The very small 
difference to the complete result multiplied with a factor 10 is plotted with dotted lines. 



the constraint y < 1, with y := p ■ q/p ■ I = Q 2 /(2xbjP • I)- The plot of da/dQ 2 is based 
on the xbj- dependent cross section integrated between x B j = 0.19 and xbj = 0.69 taking 
into account the same kinematical constraints. The dashed lines are obtained neglecting 
the contribution of the strange quark GPD and the polarized gluon GPD proving the 
dominance of the contribution of F^r, £,£). The negligible small contribution of the 
polarized gluon GPD multiplied with a factor 100 is plotted with dotted lines. 

To illustrate the dependence of the cross section on the shape of the D s distribution 
amplitude and the GPD's we show in Fig. |3| the results obtained by modeling <J> D - using 
its asymptotic form Q D -(z) = 6f£, 3 z(l — z) instead of the parameterization of Eq. @ 
(dashed lines) and alternatively by replacing the models for the GPD's F(t, £, 0) by their 
forward limit F(t, £, 0) (dotted lines). The latter choice corresponds to the approximation 
of the GPD's for t = by usual parton distributions. For comparison also the original 
result is plotted again in Fig. ||. 

In Fig. ^ we finally show the result that is obtained if all corrections of the order 
m 2 Da /Q 2 in the amplitude are neglected. The very small difference to the complete result 
multiplied with a factor 10 is plotted with dotted lines. The smallness of these corrections 
justifies the approach based on the factorization theorem, which is valid only up to terms 
of the order fi 2 /Q 2 with fi being a typical mass scale of the process. 

Integrating all variables Q 2 , xbj, and t over the kinematical region specified above 
gives a value for the total cross section of a = 2.45 x 10~ 14 GeV -2 = 9.5 x 1CT 6 pb. Even 
given this small value the huge integrated luminosities of / Ldt > 10 45 cm -2 = 10 9 pb -1 
available at a neutrino factory || would lead to a sizable number of events of the order of 
magnitude 10 4 . Larger total cross sections can be obtained for higher neutrino energies 
because of the larger available kinematical region. 

Uncertainties of the rough estimate presented here result from the lack of knowledge 
about the exact form of the meson distribution amplitude as illustrated in Fig. |3|. Also 
the predictions for the value of fo s differ by about 20%. Up to now the experimental 



uncertainty for fu B is even larger |T7|]. It is worth noting, however, that experiments of 
the kind discussed in the present article can provide much more precise information on 
the D s decay constant independently of the exact cross section for the process ([!]). As 
discussed in JL8J the relatively high production rate of D s mesons allows to determine fo s 



by measuring the D s branching ratios and its total width. 

From the experimental point of view a difficulty arises from the fact that all three 
particles in the final channel need to be detected because the neutrino beam energy is 
not sharp and therefore an exclusive event can not be identified by reconstructing the 
momentum of an undetected particle via a missing momentum analysis. 

It has been shown that the high intensity neutrino beams available at neutrino factories 
allow to study hard exclusive meson production also for weak interaction induced reactions 
opening a new possibility to study the nucleon structure and allowing to better test the 
theory of these exclusive processes. The production of charmed strange mesons proves to 
be of particular interest due to its high sensitivity to the gluon GPD's. 
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